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Abstract. We consider a Fabry-Perot cavity made by two moving mirrors and 
driven by an intense classical laser field. We show that stationary entanglement 
between two vibrational modes of the mirrors, with efli'ective mass of the order 
of micrograms, can be generated by means of radiation pressure. The resulting 
entanglement is however quite fragile with respect to temperature. 
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1. Introduction 

Quantum entanglement is a physical phenomenon in which the quantum states of 
two or more systems can only be described with reference to each other. It is now 
intensively studied not just because of its critical role in setting the boundary between 
classical and quantum world, but also because it is an important physical resource 
that allows performing communication and computation tasks with an efficiency 
which is not achievable classically [1]. In particular, both from a conceptual and 
a practical point of view, it is important to investigate under which conditions 
entanglement between macroscopic objects, each containing a large number of the 
constituents, can arise. Entanglement between two atomic ensembles has been 
successfully demonstrated in Ref. [2J by sending pulses of coherent light through 
two atomic vapor cells. More recently Ref. has shown spectroscopic evidence 
for the creation of entangled macroscopic quantum states in two current-biased 
Josephson-j unction qubits coupled by a capacitor. The interest has been also extended 
to micro- and nano- mechanical oscillators, which have been shown to be highly 
controllable [4] and represent natural candidates for quantum limited measurements 
and for testing decoherence theories [5]. Recent proposals suggested to entangle a 
nano- mechanical oscillator with a Cooper-pair box [6], arrays of nano-mechanical 
oscillators [7] , two mirrors of an optical ring cavity [5] , or two mirrors of two different 
cavities illuminated with entangled light beams [9]. These two latter proposals 
employed the optomechanical coupling provided by radiation pressure, which has 
been demonstrated to provide a useful tool to manipulate the quantum state of light 

[iniiiiiiiaiisiiiiiiiiiiiii. 

Here we study the simplest scheme in which one can test the entangling 
capabilities of radiation pressure, that is, a linear Fabry-Perot cavity with two 
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vibrating mirrors (see Fig. [T]). This system corresponds to a simplified version of 
the system of Ref. [T7], where a double-cavity set-up formed by a hnear Fabry-Perot 
cavity and a "folded" ring cavity is considered. Similarly to what has been done in 
Ref. [T7j, we determine here the exact steady state of the system and show that 
if the cavity is appropriately detuned, one can generate stationary entanglement 
between macroscopic oscillators (effective mass 100 ng). As it will be discussed 
below, the main advantages of the present scheme with respect to that of Ref. [17] 
are its simplicity and the fact that steady-state entanglement is achievable even 
with purely classical driving light, while Ref. [T7] considered the limiting case of 
large mechanical frequencies where entanglement can be generated only by injecting 
nonclassical squeezed light into the two cavities. 

The paper is organized as follows. In Section II we describe the dynamics of the 
system in terms of quantum Langcvin equations. In Section III we solve the dynamics 
and derive the correlation matrix of the steady state of the system. In Section IV we 
quantify the mechanical entanglement in terms of the logarithmic negativity, while 
in Section V we compare the present scheme with other recent proposals for the 
generation of mechanical entanglement and discuss how one can detect it. Section VI 
is for concluding remarks. 
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Figure 1. Schematic description of the system under study. The cavity (with 
mode a) is driven by a laser, and the vibrating mirrors 1 and 2 are the ones we 
want to entangle. 



2. The system 

We consider an optical Fabry-Perot cavity in which both mirrors can move under 
the effect of the radiation pressure force (see Fig. [T]). The motion of each mirror is 
described by the excitation of several degrees of freedom which have different resonant 
frequencies. However, a single frequency mode can be considered for each mirror 
when a bandpass filter in the detection scheme is used [HI and mode-mode coupling 
is negligible. Therefore we will consider a single mechanical mode for each mirror, 
modeled as an harmonic oscillator with frequency flk and effective mass Mfc, k — 1,2, 
so that the mechanical Hamiltonian of the mirrors is given by 

p2 n2 1 1 

2Mi 2M2 ^ 2 2 2 2 > \ ) 
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with [Qfej-Pj] — iTiSkj. In the adiabatic hmit in which the mirror frequencies are 
much smaller than the cavity free spectral range c/2L [L is the cavity length in the 
absence of the intracavity field) |20| , one can focus on one cavity mode only because 
photon scattering into other modes can be neglected, and one has the following total 
Hamiltonian 

H = Hrn + fiuj^a'^a + n—a^a{Qi - Q2) + ihE{e-''^''*a'< - e''^^*a), (2) 
L 

where a and ([a, a^] — 1) are the annihilation and creation operators of the cavity 
mode with frequency tOc and decay rate k, and the last two terms in Eq. ([2]) describe 
the driving laser with frequency lol and E is related to the input laser power P by 
\E\ ~ y^2PK/hujL. In general the mirror potential is also determined by the additional 
static Casimir term Vcas = — ?i7r/ [2A{Q2 — Qi + L)] [20^, which however is negligible 
for typical optical cavities with L ~ 1 cm and mirrors with effective masses in the 
/ig-ng range. 

The full dynamics of the system is described by a set of nonlinear Langevin 
equations, including the effects of vacuum radiation noise and the quantum Brownian 
noise acting on the mirrors. In the interaction picture with respect to huoLo'^ a 

a = -{K + iAo)a ~ i'^aiQi - Q2) + E + V2^a'", (3) 
Qk^Pk/Mk, A: = 1,2, (4) 
Pk = -AmlQk - "fkPk + {-ifn'^a^a + M^W^ fc = 1, 2, (5) 

where Aq — lOc ^ ojl and is the mechanical damping rate of mirror k. We have 
introduced the radiation input noise a^"{t), whose only nonzero correlation function 

is [n] 

{a'^{t)a^-^\t'))^6{t-t'), (6) 

and the Hermitian Brownian noise operators Wj(t), with zero mean value and 
possessing the following correlation functions [HI [53] 

where ks is the Boltzmann constant and T is the equilibrium temperature, assumed 
to be equal for the two mirrors. 

We are interested in the dynamics of the quantum fluctuations around the 
steady state of the system. We can rewrite each Heisenberg operator as a c-number 
steady state value plus an additional fluctuation operator with zero mean value, 
a — Us + Sa, Qk = Qi + SQk, Pk = Pj^ + 5Pk- Inserting these expressions into the 
Langevin equations of Eqs. ([3|), these latter decouple into a set of nonlinear algebraic 
equations for the steady state values and a set of quantum Langevin equations for 
the fluctuation operators. The steady state values are given by P^ — 0, {k — 1,2), 
QI = {—1)^ {hu!c/Mki^lL)\as\'^ , as = E / (n + iA), where the latter equation is in fact a 
nonlinear equation determining the stationary intracavity field amplitude as , because 
the effective cavity detuning A, including radiation pressure effects, is given by 

A = A„ + (O; - (31) (8) 



/ Hlu \ 

coth — 

\2kBT 



(7) 
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The exact quantum Langevin equations for the fluctuations are 

Sd = -(k + iA)6a - (os + Sa) {6Qi - 6Q2) + V2Ka™, (10) 
dQk=SPk/Mk fc = l,2, (11) 
SPk = -AmlSQk - ikSPk + i-^f^j; {^sSa^ + «^^«) 
+ 5aHa + MkWk k = 1,2. (12) 

From a physical point of view the strong driving regime is the most relevant one. In this 
regime, the intracavity amplitude is very large, \as\ ^ 1, and, as shown by Eqs. (fTO|l 
and ()12p . one has a large effective optomechanical coupling constant agUJc/ L between 
the fleld quadrature fluctuations and the oscillator. When |as| ^ 1, one can safely 
neglect the cavity field fluctuation operator 6a with respect to as in Eqs. PH)) and 
(|12p and consider linearized Langevin equations. Notice that this amounts to linearize 
only with respect to the cavity mode and not with respect to the mechanical oscillator, 
whose operators appear linearly in the dynamical equations from the beginning and 
therefore are not approximated in the linearized treatment. 

It is evident that the cavity mode is coupled only to the relative motion of the 
two mirrors and it is therefore convenient to rewrite the above equations in terms of 
the fluctuations of the relative and center-of-mass coordinates, i.e., 

SQcrrr - Tl^SQi + ^5Q2, SP„n - 5P, + 5P2, (13) 

where Mt — Mi + M2 and /i = M1M2/MT are the total and reduced mass of the two 
oscillators. The linearized Langevin equations for these coordinates are 

Sd = - {K + iA)Sa + i^aJQr + V2i^a''', (15) 

SQr =SPr/fl, (16) 
SPr = — fifl^SQr — JrSPr — fJ- {^2 ^ ^l) SQcm 

- (72 - 71) SPcm + (a*Ja + a,(5at) + ^Wr, (17) 

SQcm = 6Pcm/MT, (18) 
SPcm = - Mx^l^SQcm - IcmSPcm 

- M {nl - nl) 6Qr - (72 - 71) SPr + MtW,^, (19) 

where we have defined the center-of-mass frequency 51^^ — (Mi^lf + M2^'^ /Mtj 
damping rate = (M171 -I- M272) /Mt, and Brownian noise Wcm ~ 

{MiWi + M2W2) /Mt, and also the relative motion frequency ^1 = {M2^l + Mi^l) /Mt, 
damping rate 7^ = (A/27i + ^^172) /Mt and Brownian noise Wr = W2 — Wi. Thanks 
to these definitions, the center-of-mass and relative motion Brownian noise possess 
correlation functions analogous to those of Eq. ((T]) , with the corresponding damping 
rate and mass. The two noises are however correlated in general, because 

{W,,n{t)Wr{t')) = {Wr{t)W,^{t')) 

flLJ 



Mt ./_eo 27r 



coth 



2kBT 



(20) 
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The above equations show that, even though the cavity mode directly interacts only 
with the relative motion, the three modes are all coupled because of the center-of- 
mass-relative-motion coupling, which is present whenever Qi 7^ 5^2 or 71 7^ 72- 



2.1. Equal frequencies and damping rates 

The dynamics considerably simplify when Cli = Q2 = ^ and 71 = 72 = 7- In fact, 
in such a case ficm = = ^ and 7cm = 7r = 7 and the center-of-mass motion 
fully decouples from the cavity mode and the relative motion, even if the masses are 
different. The center-of-mass becomes an isolated quantum oscillator with mass Mt 
and subject to quantum Brownian noise, i.e., 

6Qc^ = SPc^/Mt, (21) 

SPcm = ~MTnHQ,„, - -fdPc,n + MrWcn,, (22) 

while the relative position of the two mirrors and the linearized fluctuations of the 
cavity mode form a system of two interacting modes described by the following linear 
Langevin equations 

5Qr = SPr/ti, (23) 

SPr = -^in^SQr - -fSPr + + flWr, (24) 



X = -kX + AY + V2nX'", (25) 

Y = -kY -AX+ ^'^/"' SQr + V2^Y"\ (26) 

where we have chosen the phase reference of the cavity field so that Ug is real, we 
have defined the cavity field quadratures X = (Sa + Sa^ )/ V2 and r = ((5a - (5at)/i V2, 
and the corresponding Hermitian input noise operators X*" = (o™ -I- a"^'^)/V2 and 
F*" = (a™ — a^"-'^)/i^/2. Notice that Eqs. (|23p coincide with the linearized equations 
of a Fabry-Perot cavity with only one movable mirror with mass fi. 

It is convenient to switch to dimensionless dynamical variables for the mechanical 
oscillators. If we define 



^Qk - \ TTTlIk: ^Pk = V^Mk^pk, fc = 1, 2, (27) 




^9cm, SPcn = ^/flMT^Pam, (28) 



SQr = W —^Qr, SPr — \/h^D,pr, (29) 



such that [qj,Pk] ~ iSjk, either for j, k = 1,2 and for j, k = r, cm, definitions p^ - (fT4| 
become 

Qcm = riqi + r2q2, Pcm = ''iPl + ^2^2, (30) 
qr =riq2-r2qi, Pr ^ rip2 - r2Pi, (31) 



where = a/ Mk/Mx, k = 1,2. The quantum Langevin equations become in terms 
of these dimensionless continuous variables 

Qcm — ^Pcm: (32) 
Pcm — 

-nq 
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Pr = Ctqr - ^Pr + GX + ^r, 

X = -kX + AY + V2i^X'"-, 
Y = -kY - ax + Gqr + \/2^F'", 
where we have defined the effective optomechanical coupling constant 



G 




(34) 
(35) 
(36) 
(37) 



(38) 



which, being proportional to the square root of the input power, can be made quite 
large, and the zero- mean scaled Brownian noise operators ^r{t) and Ccm(O) with 
correlation functions 



coth 



\2kBT) 



1 



(39) 



where j, k =r,cm. 



3. Stationary correlation matrix of the two mirrors 



When the three-mode system is stable, it reaches a unique steady state, independently 

from the initial condition. Since the quantum noises ^cm, ^r, ^™ and Y™ arc 
zero-mean quantum Gaussian noises and the dynamics is linearized, the quantum 
steady state for the fluctuations is a zero-mean Gaussian state, fully characterized by 
its 6 X 6 correlation matrix (CM) V^^^ = {ui{oo)uj{oo) + Uj{oo)ui{oo)) /2, where 
U'^(oo) = {qi{oo),pi{oo),q2{oo),p2{oo),X{oo),Y{oo)) is the vector of continuous 
variables (CV) fluctuation operators at the steady state {t oo). We are interested 
in the stationary reduced state of the two mirrors, which is obtained by tracing out 
the cavity mode. This state is obviously still Gaussian and fully characterized by the 
4x4 matrix Vij formed by the first four rows and columns of V^p . The general form 
of V is quite simple. First of all it is V12 = V34 = 0. In fact, since pj = qj/fl, j = 1, 2, 
it is 

= i<l^i^)M^)+Mo-)l^i^)) = 1 nm l{ql(t)) = 0, (40) 
2 2S2 t-too at 

and the same happens for V34. Moreover, thanks to the decoupling between center- 

of-mass and relative motion it is V14 = V23 = 0, because 

(gi (00)^2(00) -1-^2(00)91(00)) 



rir2 



Pcm 

(oo)gcm(oo)) 



Pr(oo)gr(oo))] 

2ni^-^ dt^^''-m-{qm=0, 



- {qr {00)pr{00) 

rir2 ,. d 
lim 



and the same happens for V23. The final form of V is 

V = 





Vi3 





V22 







V33 





V24 



144 / 



(41) 



(42) 
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where 

Vu = rl{ql^),t + rl{ql)su V22 = rfipU^t + rjipD.t, (43) 
^33 = rl(ql^),t + rl{ql)st. V44 = rl{pljst + rl{pl)st. (44) 
Vi3^rir2{{qUst- {qr)st) , ^ r^r2 {{plJ st - {pD st) . (45) 

that is, it depends upon the mass ratios rk and the four stationary variances 

(u^)st EE hmt^oo(w^(t)), U = qcm,Pcni,qr,Pr- 



3.1. Calculation oj the stationary variances 



The center-of-mass and relative motion stationary variances can be obtained by 
solving Eqs. ([32]) and considering the limit i — > cx). Defining the six-dimensional 
vector of variables v'^ [t) = {qcm{t) , Pcm{t) , qr{t),Pr{t), X{t),Y{t)), the vector of noises 



(46) 



n^(t) = {0,^cm{t),0,^rit),V2KX"'{t),V2KV'{t)) and the matrix 



A 



( 





n 
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Eqs. (|32p can be rewritten in compact form as v(t) — Av{t) 



v{t) = M{t)v{0) 



dsM{s)n{t- s), 



n(t), whose solution is 
(47) 



where M{t) = exp{At}. The system is stable and reaches its steady state when all the 
eigenvalues of A have negative real parts so that M(oo) = 0. The stability conditions 
can be derived by applying the Routh-Hurwitz criterion 1261, yielding the following 
two nontrivial conditions on the system parameters. 



si = 27K [A"* + A^{-/^ + 27K + 2^2 - 



(48) 



+ (7K + + fl^f] + nG^A{-f + 2Kf > 0, 

S2^n (A^ + K^) - G^A > 0, (49) 

which will be considered to be satisfied from now on. If we consider the variables 
Vj{t)., we can construct the stationary correlation matrix 

{vi{oo)vj{oo) + Vj{oo)vi{oo)) 



Cij — 



which is the quantity of interest because Cn — {q'^^)st, C22 — {p\ 
and C44 = {Pr)st- When the system is stable, using Eq. ([47]) one gets 



2 \ 

cm/ St J 



C33 — 



poo 1^00 

ds'A'hk{s)Mji{s')<^ki{s 

k,i 



(50) 



(51) 



where $fci(s — s') = {{n]^{s)ni[s') + ni(s')nk{s))) /2 is the matrix of the stationary 
noise correlation functions. Due to Eq. (j39p . the mirror Brownian noises are not delta- 
correlated and therefore do not describe in general a Markovian process. However, as 
we shall see, mechanical entanglement is achievable only using oscillators with a very 
good mechanical quality factor Q = ^l/^. In this weak damping limit, 7^0, the 
quantum Brownian noises ^r{t) and ^cm{i) become delta-correlated, [27j 



{m^k(t') + Ut')^M /2 - 5,kl i2n + 1) 5{t - t'), 



(52) 
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where n = {ejip{hn/ kgT} — 1) ^, is the mean thermal excitation number, and one 
recovers a Markovian process. Using the definitions ol Xj„ and Yi^ and Eq. we 
finally get $fe;(s - s') = DkiS{s - s'), where 

/ \ 

7(2n+l) 

K ■ 

\ Ok/ 

As a consequence, Eq. ([STj) becomes 



(53) 



C = 



dsM{s)DM{s) 



(54) 



which, when the stability conditions are satisfied so that Af(oo) = 0, is equivalent to 
the following equation for the CM, 

AC + CA^ = ~D. (55) 

Eq. (|55p is a linear equation for C and it can be straightforwardly solved. The center- 
of-mass is decoupled from the other two modes and Eq. fSS]) trivially gives 



1 



Cll = C22 (gcm>st = (Pcm>st = o + 



(56) 



The relative motion is instead coupled with the cavity mode and consequently the 
final expression of the stationary variances are much more involved. One has 



C33 
C44 



{ql)st 



dpH, 



(57) 
(58) 



where 



-^G^K {A2(7 + k) + K{-fK + + n^) - Af7(7 + 2k)} , (59) 
dp = l-[si]-'^2G'^KnA{j + 2K), (60) 
bq = [2siS2r^G^ |2k (A^ + K^) I [a2 + (7 + k)2] {k^ + 7A) (61) 

+ f^^ (7 + k) {n - 2A) } AG^r^^ ^ 2k) [A7 - K (f7 - 2A)] } , 
dq^l + [siSal^^AG^ [si - 27Kr22 (f]2 27K + 4k2) 



4K^f^^ A^ 



(62) 



4. Conditions for stationary entanglement 

Simon's separability PPT (positive partial transpose) criterion is necessary and 
sufficient for bipartite Gaussian CV states [2^. It assumes a particularly simple 
form for the CM of the two mirrors of Eq. (|42p . In fact, after some algebra, one 
gets the following necessary and sufficient condition for the presence of mechanical 
entanglement between the two mirrors in the stationary state. 



\Qcm)st {Pcm)st 



{pl)st{qcm)st ~ 
(ll) st(j>l) St 



< 1- 



1 



(63) 
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where we have defined = 4rjr| — A^/M^- For very different masses ?/ — > and 
the right hand side of Eq. ([551) tends to —oo, i.e., the criterion is never satisfied and 
the mirrors are never entangled. It is evident therefore that stationary entanglement 
is better achieved for equal mirrors, i.e., = 1, when the right hand side of Eq. (j63p 
is equal to zero and the necessary and sufRcient entanglement condition becomes 
equivalent to a "product" of sufRcient criteria analogous to those derived in [H [25] , 
that is, {qf)st{plm)st < 1/4 or (Pr) stiQcm) st < 1/4. Since the center-of-mass of the two 
mirror is unaffected by the optomechanical coupling (see Eq. l[56[) ). this means that 
the two mirror vibrational modes are entangled if and only if their relative motion is 
sufficiently squeezed, i.e., 

{l'r)s. or {pD., < ^^^1^. (64) 

This equation provides an intuitive picture of how the entanglement between the two 
mirrors is generated by the radiation pressure of the light bouncing between them. 
If the cavity is strongly driven, the radiation pressure coupling becomes very large 
and the fluctuations of the mirror relative motion can be significantly squeezed. If 
such a squeezing is large enough to overcome even the thermal noise acting on the 
center-of-mass, Eq. ([64[) guarantees that the two mirrors are entangled. Eq. ([64[) 
also points out the main limit of the proposed scheme: the mirrors center-of-mass is 
not affected by radiation pressure and cannot be squeezed. This suggests that the 
generated entanglement is not robust against temperature because satisfying Eq. ([64p 
becomes prohibitive at large n. 

One can quantify the stationary mechanical entanglement by considering the 
logarithmic negativity Ej^ 28J, which in the CV case E_^ can be defined as [29] 

E^r = max[0,-ln2i^"], (65) 

where v~ is given by 



2-1/2 



S(y)- (S(l/)2-4detl/)'^' (66) 



1/2 



with S(y) = det Ni + det N2 — 2 det N12 and we have used the 2x2 block form of the 
CM 

Therefore, a Gaussian state is entangled if and only if i/^ < 1/2, which is equivalent 
to Simon's necessary and sufRcient entanglement criterion for Gaussian states [24] of 
Eq. ([55]) . and which can be written as 4dety < S(F) — 1/4. In the case of the 
stationary matrix V of Eq. ([1^ , one has 

detV= {ql)st{pl)st{qlJst{plJst (68) 
= {l-ri) [{ql)st{pl)st + {qlJst{plJst\ 

+ ^ [{qD st{pl,n) St + {ql,ri) st{pl) St] ■ (69) 

Therefore, in the most convenient condition for entanglement, i.e., identical mirrors 
<^ ?7 = 1, one has Y.{V) ^ {qf)st{pUst + {qUst{pDst, yielding 

V' = V('?r>st(p?m)st, \/ {qlri)st{pl)st'^ , (70) 
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so that in this case of equal masses, the logarithmic negativity assumes the particularly 
simple form 

E'aa = max|o, -In 2^/ {ql)st{pl^)st , - In 2yJ {ql^^} st{pl) st } -(71) 

Using Eqs. and ([7T|) one has stationary entanglement if one of the two 

following conditions is satisfied 

n 

b, + d,n<-^^, (72) 

bp + dpn<- ^_^^ , (73) 

which, as expected, are better satisfied in the zero temperature limit, n — s- 0, since 
dq,dp > whenever the stability conditions are satisfied (otherwise one could have 
negative variances at large enough temperatures). 

These two equations lead us to the main result of the paper, i.e., it is possible 
to realize an entangled stationary state of two macroscopic movable mirrors of 
a classically driven Fabry-Perot cavity. However, such a stationary mechanical 
entanglement turns out to be fragile with respect to temperature, as it can be easily 
grasped from Eqs. (|72p - (|75|) . This is illustrated in Figs.[2]131 where we have considered 
a parameter region very close to that of recently performed experiments employing 
optical Fabry-Perot cavities with at least one micromechanical mirror [TTJ [121 ESI Ell ■ 
Figs. |2]|3l refer to the case of an optical cavity of length L — 1 cm, finesse T = 1.9 x 10^, 
so that K ~ 5 X 10^ s~^, driven by a laser with wavelength 1064 nm and power P = 50 
mW. The two identical mechanical oscillators have angular frequency il/2TT — 10 
MHz, damping rate 7 = 3 x 10^ s^^, and mass m — 100 ng. Fig. [51 refers to the zero 
temperature limit and shows that stationary entanglement is present only within a 
small interval of values of A around A ~ Aopt where 

This value is essentially the optimal value for the detuning for achieving entanglement. 
This can be understood from the expression of {p'^)st- In fact, at zero temperature 
entanglement is obtained when bp < (see Eq. ((73l) ). which is satisfied when the 
numerator of Eq. (|59p is negative, since si > due to stability. This condition 
is obtained by considering the minimum of the second order polynomial in A in the 
numerator and by imposing that it is negative. The minimum value of this polynomial 
is obtained just at A = Aopt and it is negative when 

7f7>2K(7 + K). (75) 

Therefore A = Aopt and Eq. ([TSl are sufficient conditions for achieving entanglement 
at zero temperature. This parameter regime is the optimal for entanglement because 
when A ~ Aopt — si is also close to its minimum value, implying therefore a 
large negative value of bp (see Eq. ([59l)) and also a value of dp very close to zero 
(see Eq. (j60p ). which means an improved robustness of entanglement with respect 
to temperature. In Fig. |3l we study the resistance to thermal effects by plotting Ej^ 
evaluated at the optimal detuning, i.e., corresponding to the maximum of Fig.j^l versus 
temperature. We see that this entanglement vanishes for T > 100 /iK. This behavior is 
valid in general, even in parameter regions different from that of Figs. [31 whenever 
one finds a regime with a nonzero stationary entanglement, this entanglement quickly 
tends to zero for increasing temperatures. As discussed above (see below Eq. (|64|) this 
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is due to the fact that in this simple Fabry-Perot cavity system, the mirror center- 
of-mass is unaffected by the radiation pressure of the cavity mode and remains at 
thermal equilibrium. One could achieve a larger and more robust entanglement by 
adopting the double-cavity setup considered in |17j . where the optical mode of the 
second, "folded" cavity couples just to the center-of-mass of the mirrors of interest, 
which is then also squeezed, independently from the relative motion. In this latter 
scheme therefore robustness against temperature is achieved at the price of a much 
more involved experimental setup. 

Eq. (|7ip shows that mechanical entanglement at zero temperature could be 
realized as well when ((?^)st < 1/2. However, it is possible to see through numerical 
calculations that this condition is much more difficult to realize with respect to 
{Pr)st < 1/2. This fact is not easily seen from the analytical expressions of bg and dq 
(Eqs. (f6T |) - (f62|) ). which are more difficult to analyze with respect to those of bp and 
dp (Eqs. dSSD-dini))- 




Figure 2. Logarithmic negativity Ej\f as a function of the normaUzed detuning 
A/f! in the case of an optical cavity of length L = 1 cm, finesse = 1.9 X 10^, 
driven by a laser with wavelength 1064 nm and power P = 50 mW. The identical 
movable mirrors have a frequency Q/2tt = 10 MHz, damping rate 7 = 3 X 10* 
s~^, mass m = 100 ng and their temperature is T = 0. 



5. Comparison with other proposals and experimental detection of the 
entanglement 

It is interesting to compare the present proposal with other recent schemes for 
entangling two micro- mechanical mirrors, especially with Refs. |S1 [HI HZl US] ) which 
are all based on the optomechanical coupling provided by the radiation pressure. 
Refs. [U [9l [17] considered the steady state of different systems of driven cavities: 
Ref. [8] focused on two mirrors of a ring cavity and considered the situation in 
the frequency domain; Ref. [9] assumed to drive two independent linear cavities 
with two-mode squeezed light and stationary mechanical entanglement is achieved 
by transferring the entanglement of the two driving beams to the end-mirrors of the 
two cavities. 
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Figure 3. Logarithmic negativity Ej\f evaluated at A = Aopt — O.SQ (see 
Eg. I74I I) as a function of temperature. The other parameter values are the same 
as in Fig. |2] 



We have already partially compared the present scheme with that of Ref. |17j . 
to which is strongly related. In fact, the double-cavity scheme of Ref. [17] coincides 
with the single Fabry-Perot cavity scheme considered here when the "folded" cavity of 
Ref. [17] is not driven. The additional folded cavity couples to the center-of-mass of the 
two vibrational modes and if it is appropriately driven by squeezed light, it is able to 
transfer this squeezing to the center-of-mass. This has the advantage of increasing the 
entanglement and making it more robust against temperature (see Eq. (ITTI) '). but this 
is obtained at the price of a more involved apparatus, requiring the preparation of an 
additional ring cavity and the use of nonclassical driving. Moreover, Ref. [17j evaluated 
the stationary state of the two mechanical modes approximately, by considering the 
resonant case A = and solving the dynamics of the system only in the limit when 
fl is much larger than the other parameters, G, k, so that fast terms rotating at 
frequency can be neglected in the equations of motion. In this limit, Ref. |17j finds 
that the steady state of the two mirrors is entangled only if the input field is squeezed, 
while is never entangled for a classical coherent input. Here we determine the steady 
state of the system exactly in the Markovian limit of weak mechanical damping and 
we find that when fast terms rotating at frequency f2 cannot be neglected, one can 
entangle the mirror even using classical driving. As expected, it is possible to check 
that the present exact solution reproduces the results of Ref. [17j in the same limiting 
conditions (no input squeezing, large mechanical frequency, and no folded cavity). In 
fact, if we consider A = 3> G, k, 7 in Eqs. (|59ll - ([62|) . one gets 



0, 



7 (G^ + 2jK + 4k2) 
(7 + 2k) (G2 + 27K)' 



so that 



{qr)st = {pI)s 



1 7 (G^ -f 27K + 4k2) 

— -1-77 1 — 

2 (7 + 2k)(G2 + 27k) " 2' 



> 



(76) 
(77) 

(78) 



coinciding with Eq. (23) of Ref. [17] in the case of no input squeezing, and implying 
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absence of mechanical entanglement. Therefore we see that the "resonance" condition 
A = J7 is very close to the optimal condition for generating mechanical entanglement, 
and that if one leaves the regime of vary large mechanical frequencies ^ G, k, 7, 
one can achieve stationary mechanical entanglement even without input squeezing. In 
fact, the parameter regime considered in Figs. [2][3] corresponds to ~ G 3> k — 7. 

Another recent proposal employing radiation pressure effects for entangling two 
vibrating micro- mirrors is Ref. 18\, where the radiation pressure of an intense laser 
field first generates optomechanical entanglement between a mirror vibrational mode 
and an optical sideband. Such an entanglement is then swapped to two separated 
micro-mechanical oscillators via homodyne measurements on the optical modes, 
representing Bell measurements in this continuous variable setting. In this latter 
proposal, macroscopic mechanical entanglement is generated when the homodyne 
measurement is performed and it is therefore a transient phenomenon, with a lifetime 
which is severely limited by the mirror thermal reservoir [TH] . In the present scheme, on 
the contrary, mechanical entanglement has an infinite lifetime because it is generated 
at the steady state, and therefore its experimental detection becomes much easier. 

We also notice that the system studied here is similar to the one considered in 
Ref. [15j , where a Fabry-Perot cavity with only one vibrating mirror is considered. In 
Ref. [H] a quantum Langevin treatment analogous to the one adopted here is used to 
quantify the amount of bipartite entanglement between the vibrational mode of the 
mirror and the intracavity field at the steady state of the system. 

We finally discuss the experimental detection of the generated mechanical 
entanglement. The measurement of Ej^ at the steady state is quite involved because 
one has to measure all the ten independent entries of the steady state correlation 
matrix V . This has been recently experimentally realized (see Ref. |30[ and references 
therein) for the case of two entangled optical modes at the output of a parametric 
oscillator. Instead, one does not have direct access to the vibrational modes and 
therefore it is not clear how to measure them. However Ref. [H] showed that, apart 
from additional detection shot noise, the motional state of the mirror can be read from 
the output of an adjacent Fabry-Perot cavity, formed by the mirror to be detected 
and another "fixed" (i.e. with large mass) mirror. In fact, it is possible to adjust the 
parameters of this second cavity so that both position and momentum of the mirror 
can be experimentally determined by homodyning the cavity output light |15j . In 
particular, if the readout cavity is driven by a much weaker laser so that its back- 
action on the mechanical mode can be neglected, its detuning A2 is chosen to be 
equal to the mechanical frequency 17, and its bandwidth K2 is large enough so that 
the cavity mode adiabatically follows the mirror dynamics, the output of the readout 
cavity a™* is given by 

a°"* = i-^fo + a^, (79) 

where h is the annihilation operator of the vibrational mode, G2 is the effective 
optomechanical coupling of the readout cavity (see Eq. (P5|) ). and a™ is the input 
noise entering the readout cavity. Therefore using a readout cavity for each mirror, 
changing the phases of the two local oscillators and measuring the correlations between 
the two readout cavity output one can then detect all the entries of the correlation 
matrix V and from them numerically extract the logarithmic negativity £V by means 
of Eqs. and 
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6. Conclusions 

We have considered a system formed by a linear cavity with two vibrating mirrors, 
driven by an intense classical light field. The two mirror vibrational modes interact 
thanks to the radiation pressure of the light bouncing between them. We have 
determined the steady state of the system and we have seen that, in the case of identical 
mechanical oscillators, the two vibrational modes become entangled if the cavity 
detuning is close to the mechanical frequency. The resulting mechanical entanglement 
is however quite fragile with respect to temperature and this suggests that, in order 
to generate macroscopic mechanical entanglement which is more robust with respect 
to thermal effects, it is convenient to drive the cavity with nonclassical light (see e.g., 

m-) 
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